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Abstract 

Theory of scattering by many smah bodies is developed under vari- 
ous assumptions concerning the ratio ^, where a is the characteristic 
dimension of a small body and d is the distance between neighbor- 
ing bodies d = 0{a'^^), < ki < 1. On the boundary Sm of every 
small body an impedance- type condition is assumed un = CmU on Sm, 
l<m<M,Crn = h , < K, hfn are constants independent of a. 
The behavior of the field in the region in which M = M{a) ^ 1 small 
particles are embedded is studied as a — > and m(a) — > oo. Formulas 
for the refraction coefficient of the limiting medium are derived under 
the assumptions: a) ki = (2 — k)/3, < k < 1, and b) ki = 1/3, 

K> 1. 
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1 Introduction 

The theory of wave scattering by small bodies was originated by Rayleigh 
in 1871 [1]. In ^ this theory was developed for small bodies of arbitrary 
shapes, analytic formulas for the S'-matrix for acoustic and electromagnetic 
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(EM) wave scattering by small bodies of arbitrary shapes have been derived. 
These formulas allow one to calculate the S*- matrix with any desired accuracy. 
Analytic formulas for the electric and magnetic polarizability sensors have 
been derived for bodies of arbitrary shapes [2], |3]. In [3J - yjj a theory 
of wave scattering by many small bodies embedded in a bounded domain 
filled in by a material with known properties was developed. It was assumed 
in P and (TU] that the characteristic size of the small particles (bodies) is 
a, that the distance d between two neighboring particles is is of the order 
d = 0{a^^^), that the total number of the embedded particles M = 0{^), 
and that the boundary condition on the boundary Sm of m-th particle 
is of impedance type: 

un = Cmu on Sm, 1 < m < M, (1) 



where is the unit normal to 5*^ directed out of Dm, and Cm = ^ , where 
hm, Im hm < 0, is a constant independent of a. 

The waves in the original material are described by the equation 

Louo := [V^ + kVQ{x)]uo = in (2) 

where 

nl{x) = 1 inD' = R^\D, (3) 

D is a bounded domain, and UqIx) is continuous in D (or piecewise- 
continuous with a finite number of discontinuities, which are smooth sur- 
faces), ImnQ > 0. The scattering solution to ([2]) satisfies the radiation con- 
dition 

Mo = e'*'"'^ + fo, (4) 
^-^ — ikvo = o(-), r := Ixl oo. (5) 

If small particles are embedded in then the scattering problem consists of 
finding the solution to the following problem: 

M 

LoUM = in y Dm, (6) 

m=l 

= CmUM on Sm, l<m< M, (7) 
um = Uo + vm, (8) 
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where Uq solves problem ([2]), (jl]), ([2]) and vm satisfies the radiation condition 
similar to 

It is proved in fTO] that problem ([6]) - ([8]) has a unique solution and this 
solution is of the form 

M „ 

um = Uo{x) + / G{x,t) am{t)dt, (9) 

™ — 1 ^ Sm 



m=l ' 



where G{x, y) is the Green function of the operator Lq for M = 0, i.e., in the 
absence of small particles: 

LoG{x,y) = -5{x-y) inW", (10) 



G satisfies the radiation condition and am solves the equation 

UcN - CmUe + —^^^--^ — - CmTm(Tm = OU 5*^. (H) 

Here is the effective field acting on the m-th particle: 

Ue{x) := uf'\x) := UMix) - / G{x,t) a{t)dt, x G M^ (12) 



dGis t) 

Am(Tm-='2 I Om{t)dt, TmOm ■= / G{s,t)crm{t)dt. (13) 



It was proved in [TOj that 

G{x, y) = 1 [1 + 0{\x- y\)\ , \x-y\^0, (14) 

47r|x — y\ 

and one can differentiate formula (|T^ . 

The following result is also proved in [TU]. Assume that is a ball 
of radius a centered at a point x^- Let h{x) be an arbitrary continuous 
function in D, lmh{x) < 0, Ap C D be any subdomain of D, and J\f{Ap) be 
the number of particles in Ap. Assume that 

J\f(Ap) = -[ N(x)dx\l + o(l)], a^O, (15) 
a J A, 

where N{x) > is a given continuous function in D. Let 

, , 4:7iN(x)h(x) , , 

Finally, assume that Cm := ^^^^^ . Now the result can be formulated: 



Theorem 1 ([ID])- Under the above assumptions there exists the limit 

\im\\ueix) - u{x)\\c{D) = 0- (17) 

The function u{x) solves the problem 

Lu := [V^ + A;^ - q{x)\u = m (18) 
u = Uq + v, (19) 

where Uq satisfies equations (jl]), the function v satisfies the radiation 
condition siilar to and 

q{x) := qo{x) + p{x), (x) := 1 - k^"^ q{x) , (20) 

where p{x) is defined in (fT6|) . 

qo{x) := e - enlix), (21) 

and n^x) is the coefficient in ([2]). 

The aim of this paper is to investigate the behavior of Ue{x) when the as- 
sumptions Cm = ^^^"^ , d = 0(a-^/^), M = O(^) are replaced by the following 
more general assumptions: 

C„ = ^, rf = 0(a-), M = 0(J^), (22) 

where n > —1 and < ki < 1 are parameters. 

If Ki = 1, then the distance between neighboring particles is of the order 
of the size of a small particle. This is a special case which is not covered by 
a rigorous theory. However, if ki is close to 1, then practically the distance 
between neighboring particles is very close to the order a of the size of a 
small particle. 

In [TU] the theory was developed in detail in the case /t = 1, ki = |. 
The questions we are interested in this paper are: 

1) For what ranges of k, and ni the limit u{x) of Ue{x), as a — 0, does 
exist? 

2) What is the equation which this limit u{x) solves? 
The answers we give are: 



4 



1) If K < 1 and Ki = , then the hmit (JTTll exists, 



(1+0(1)), = -ATTh{x^)a'^ "UeiXm) 



and the hmiting function u{x) solves the following equation: 

u{x) = Uo{x) -471 [ G{x,y)h{y)N{y)u{y)dy. (23) 



D 



Therefore, u solves equation (|T8l) with given by (!20i) . go (a;) given by 
dnD, and 

p{x) = 47Th{x)N{x), (24) 
where A^(x) > is defined by the formula: 

Af{Ap) = -^i N{x)dx[l + 0(1)], a ^ 0. (25) 
2) If K > 1 then, as a — > 0, 

= -4'KUe{Xm)a{l+ oil)), Gm = ^-^IL. h + o{l)) , 

and the limit (ITTI) exists if /ti = | . The limiting function 'u(x) solves 
equation (l23l) with = 1 and A^(a;) defined by (!25l) . The function 
also solves equation (fT8|) with q{x) given by (120!) . go(2^) given by 
(12Ti) . and given by (1211) with = 1. 

In both cases, < k < 1 and k > 1, we have ki < 2/3. This implies that 
the total volume of the embedded particles tends to zero as a ^ 0. Indeed, 
the order of the total number of the embedded particles is 0{a~^'^^), and 
the total volume of the embedded particles is of the order 0{a^~^'^^) ^ as 
a ^ 0. 



2-K 

•II hi = 

K = —1 and Crn = h{xm)a- Moreover, one has: 



Let us make a remark about the case when ki = = 1 . In this case 



4:71 

Qm ~ [— AMe(x„) - 47lh{Xm) Ue{x^)\ ; 
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where A = is the Laplacean, and 



The quantity := \G{x,ym)Qm\ = 0[a^ '^^), as a ^ 0, and 

Jm = \ [ [Gix, t) - G{x, Xm)] <ym{t)dt\ = O (a^-^-ifl^) , a^O. 



For the relation Jm <^ Im to hold as a — > 0, it is sufficient that the relation 

holds. For this relation to hold it is sufficient to have ki < 1. 

The relation -C Im allows us to use formula (35), see below, i.e., 
approximate the exact formula (31) by an approximate formula (35) with an 
error which tends to zero as a — > 0. 

Assuming ki < 1, one has 

M , 

A-7T 

«e(a:) = Wo + ^G(x,i/(^))[— Awe(y(^))-47r/i(y(^))tz,(y(p))]aW(A,). (26) 
We have: 

aW(Ap) = ^ / N{x)dx[l + 0(1)] ^ a'-'^^ iV(y(^))|A,|, (27) 
a Jap 

where o(l) tends to zero as a ^ 0. For the limit of the sum in (l26l) to exist 
as a — > 0, it is necessary and sufficient that 3 = 3ki, i.e., ki = 1. If ki = 1, 
then the limit of Ue{x), as a — > and maxp diamAp 0, is the function 
u{x), which solves the equation 

u{x) = uo{x) + G{x, y) [— Au^y) - Anh^y) u{y)\ N{y)dy. (28) 

Applying operator Lq = + A;^ — qo{x) to (l28l) and using equation (fTOl) . one 

gets 

Lou = -[—Au-ATTh{x)u{x)]N{x). (29) 
3 
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m=l 



Thus 



4:TT 

[IH N{x)] V^M + k'^ u — qo{x)u — A'Kh{x)N{x)u{x) = 0. (30) 
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The solution u9x) to equations (18) or (30) is a locally Hf^^ function, where 
Hf^^ is the Sobolev space of twice differentiable in L^— sense functions on 
every bounded open subset of M^. This local smoothenss: u e Hf^^{R^), 
follows from known results on elliptic regularity, provided that the coefficients 
qo{x) and h{x)N{x) are in Lf^^. If these coefficients are smoother, then u is 
smoother. 

The assumption ki < 1 allows us to prove that formula (35) of Section 2 
is a good approximation of m as a — 0. 

The conclusions, obtained under the assumption ki = 1 are not proven 
to be exact in the limit a — 0. 

In Section [2] we prove the results listed in the answers. 



2 Proofs 

In the proofs we use some arguments from fTOj. 

Case 1). Consider first the case k < 1. Let us write the exact formula ([9]) 
as follows: 

M M „ 

um{x) = uo{x) + ^ G{x, Xm)Qm + ^ '^(^' t)am{t)dt, (31) 

m=l m=l'^^™- 

where Xm is the center of the ball Dm and 

Qm:= I amit)dt. (32) 



One has the following estimates (see [TU]): 

\G{x,y)\<^^, |VG(a;,2/)| <cmax(— , ^ ., ), (33) 
\x — v\ \\x — v\ [x — vV/ 



where c > stands for various constants independent of a. 

Let us estimate Qm- Eventually we want to derive sufficient condition for 
the relation 

Im ■■= \G{x,Xm)Qm\ > / [G{x,t) - G{x,Xm)](^m{t)dt := (34) 
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to hold as a — s> and \x — Xm\ ^ cl- This relation allows one to rewrite the 
exact formula fl31l) as an approximate formula: 



M 



um = uo{x) + ^ ym)Qm, \x - Xm\ > o, (35) 

m=0 

the error of which tends to zero as a ^ 0. 

To derive a formula for Qm, integrate (ITT!) over Sm and use the divergence 
theorem to get: 

-vra Aue(x„) — Ue[Xr,i)ATxa =Q„^ + — dx — -. (36) 



One has 



^ ^7^=/ rfta„(t)/" -/^ = aQ^. (37) 
' 47r s — t q q 4TT\s — t\ 



Sm ^ Si 

Here we have used the formula 

ds 



Sm- — ■ I I — c,} 47r I s i I 



a, t E Sm- (38) 



Thus, (ES]) yields: 

= ]^ J • i3yj 

If K < 1 and a — i> 0, then ( l39i) implies 



^ vra^ AMe(a:m) - 47r/i(a;m)Me(a;m)a^ 
1 + /la^-'^ 



= -47r/i(a;„)Me(x„)a'"^[l + o(l)], a ^ 0. (40) 

This is the formula for Qm which we wanted to derive. 

If a -C 1, then a formula for am can be derived as follows. The function 
Ue does not change at a small distance of order a. Therefore one can assume 
that in a neighborhood of Dm the function Ue is a constant, and one considers 
a static problem of finding am- 

Um[x) = Ue{Xm) + / "Tn IT ' ^ITT = CmMAf OU Sm- (41) 

q A7r\x — t\ oN 
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We look for the solution am = ca* ■, where c and 7 are constants. In this case 
we have: 

am{t)dt ^ 1^ _ ^ ^^^^^ ^^2) 



47r|s — t| \x — Xr.^, 

Using the impedance boundary condition (14T!) and choosing the origin at the 
point Xm, one gets 



—ca 



From the above equation one derives: 

CO? = — r . (43) 

a^'la + h{xm)a^~'^] 

If K < 1 and a — >■ 0, then equation fH3|) implies 7 = — k and c = 

-/l(Xm)Me(a;m), SO 

h{Xm}U(.{Xm) r, I /I xi n /' /l /l ^ 

c^m = 1 + 0(1), a^O. (44) 

a** 

This is the formula for am whih we wanted to derive. 

Let us find sufficient conditions for the relation (l3l|) to hold. Using esti- 
mates (l33l) and (139|) . we get 

2 — K 

G(a;,x^)Q^ = 0(^^;^), a ^ 0, |a; - a;^| > d = ©(a^^^). (45) 
Using fl33|) and (jH]) one gets: 

= ^(^ =^(«'"''""'')' k-a^-l>^- (46) 

For (IMI) to hold it is sufficient to have: 

^3-«-2«i ^ a^O. (47) 

This relation holds if ni < 1. 

Thus, let us assume that k < 1 and ki < 1, and use formulas (136!) . (139!) 
and (HO]) to get 

M 

MM(a;) = Uo{x) - ^ (^(x, Xm) 4:nh{xm)ue{xm)a^~'^, \x-Xm\ > d. (48) 



Now we want to pass to the limit a ^ in equation fl48l) . To do this, let 
us partition the domain D into a union of small cubes Ap centered at points 
y^'P'> and having no common interior points. The side of Ap is 6 ^ a. The 
number //{Ap) of small particles in Ap by formula fl25l) is: 

Ar(Ap) = ^ [ Nix)dx[l + 0(1)] = iV(y(^)) I Apl [1 + 0(1)], (49) 

where o(l) in the second equation tends to zero as diamAp tends to zero, 
and I Apl is the volume of the cube Ap. Write the sum in equation (HHI) as 

5^G(x,i/(^))47r/i(y(^))«e(l/(^))a'''^ J2 ^ 

= G(x, y^^^) 47r/.(i/(^)) w,(y(^)) ^ iV(i/(^)) I Apl (1 + 0(1)) . (50) 

a 

p 

The sum in fl50l) is a Riemannian sum for the integral 

f G{x,y)Aixh{y)N{y)u,{y)dy. (51) 

J D 

The limit of the sum in fl50|) . as a ^ 0, exists if and only if 2 — k = Ski, i.e., 
Ki = (2 - k)/3. Note that if < k < 1, then < ki < 2/3. 
In the region \x — Xm\ > one has: 

\uM{x)-Ue{x)\<o(^-^^=o{l), a^O. (52) 

The number of small particles in a unit cube is O(^) = C'(^^) if = 
0(a'^^), where d is the distance between two neighboring particles. 

We assume that the functions h{ii), Ue{y) and G{x, y) are continuous func- 
tions of so the error of replacing, for example, h{ym) by h{y^P^), where 
ym e Ap, goes to zero as diamAp 0. The function G{x,y) is not con- 
tinuous as 1/ — > X, but G{x,y) is absolutely integrable, so one may remove 
a small neighborhood of the singular point x in the integral flSTl) and the 
change of this integral will be negligible if the neighborhood is sufficiently 
small. The function h is at our disposal, and we choose it to be continuous. 
The continuity of Ue and of its limit u follows from the relation u G Hfg^{M.^). 
A more detailed argument is given in flU\ . 
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Assuming that ki = (2 — k)/3 and passing to the hmit a — in (1481) 
yields equation (12^ . Applying to this equation the operator Lq and using 
equations ([2]) and ffTOl) . one gets equation f|T8l) with q{x) given by fl20l) and 
p(a;) given by fl2^ . 

We have proved all the claims in the answer to question 1). □ 



Note that ii k < 0, then the impedance parameter ha tends to zero as 
a ^ 0. 



Case 2). Let us justify the answer to question 2). We assume now that 
K > 1. Then fl5^ implies 

Qm = -4:1TUe{Xm)a[l + o{l)], 0, (53) 

and equation fH5]) yields 

ca ' = , (54) 



so 



and 



7 = -1, C=-UeiXm), 

a^ = -^^^ [(1 + 0(1)], a^O. (55) 

Let us check when the relation ( l34l) holds, i.e., when formula (l35i) is valid, 
in other words, when formula ( l35l) yields an accurate approximation of mm, 
defined by formula (13T|) . From (!33|) and (!53|) we conclude, using the relation 
d = 0(a''i), that 

|G'(a;,xJQ^| =0(^) =O(ai''^0, a^O, \x - x^\ > d. (56) 
Furthermore, using fl55|) and fl33l) . one gets: 

J<0(^a2-i) =0(a2-2«i), a^O. (57) 
The relation holds if a^^^^^ < a^-^^S that is, if Ki < 1. 
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Let us assume that ki < 1, so that formula ([2SD is apphcable. We repeat 
the arguments given below formula fl48l) . Due to formula (1531) . now formula 
( H8|) takes the form: 

M 

Um{x) = Uo{x) — G{x, Xm) 4:1TUe{Xm)a. (58) 

m=l 

We conclude from this formula that Ue{x) tends to the limit u{x), and u 
solves the equation: 

u{x)=Uo{x)- / G{x,y)4:7iN{y)u{y)dy, (59) 
Jd 

provided that 

/^i = ^ , (60) 

and N[x) is defined by the formula fl2^ for any subdomain A G D, where 
7V(A) is the number of small particles in A. 

Applying the operator Lq to (!59|) one gets equation (ITSl) for u{x), with 
q{x) given by (l20l) . qo{x) given by (12T!) . and p{x) given by the formula 

p{x) =4:nN{x). (61) 

Since A^(a;) > 0, the function p{x) is nonnegative. 

The assumption k > 1 leads to the equation ( fT8|) with the potential g(x) 
which can vary much less than in the case k < 1, because the function h{x) 
does not enter in the definition of q{x) when k > 1. 



3 Creating materials with a desired refrac- 
tion coefficient 

If K < 1 and = (2 - /t)/3 , then equations ([T8]), ((201) and (El) hold. Thus, 
given nl{x) and n^{x), one calculates 

p{x) = k'^[nl(x) — n^(x)] := pi{x) + ip2{x). (62) 

From ([23]) and fl62|) one gets an equation for finding h{x) := /ii(x) + ih2{x) 
and A^(x) > 0: 

4:7r[hi{x) + ih2{x)]N{x) = Pi{x) + ip2{x). (63) 
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Thus 

N{x)h{x) = ^ , Nix)h2ix) = ^ . (64) 

There are many solutions {hi,h2,N} of two equations ( 16^ for the three 
unknown functions hi, /i2, N{x), /12 < 0, > 0. The condition Imn^(x) > 
imphes Imp = p2 < 0, which agrees with th inequahties h2 < 0, N > 0. One 
takes N{x) = hi{x) = h2{x) = at the points at which pi{x) = P2{x) = 0. 
At the points at which \p{x)\ > 0, one may take 

Ar(a:) = AT = const , h,{x) = ^ , h2{x) = ^ . (65) 

Let us partition D into a union of small cubes Ap, which have no common 
interior points, and which are centered at the points and embed in each 
cube Ap the number 

Af{Ap) = \^ I N{x)dx] (66) 
La Jap J 

of small balls of radius a, centered at the points Xm, where [c] stands for 
the integer nearest to c > 0. Let us put these balls at the distances 0(a~3~), 

and prepare the boundary impedances of these balls equal to ^^^^^ . 

Then the resulting material, which is obtained by embedding small par- 
ticles into D by the above recipe, will have the desired refraction coefficient 
n^(x) with an error going to zero as a — 0. 



4 Conclusions 

Wave scattering by many small particles, embedded into a material with a 
known refraction coefficient, is studied in this Letter under various assump- 
tions about the orders ki and — /t with respect to powers of a of the dis- 
tances d = 0(0"^) between the neighboring small particles and their bound- 
ary impedances ( = 0{a~'^). 

If a is the characteristic size of a spherical small particle Dm, d = 0{a'^^) is 
the distance between the neighboring particles, Cm = ^^^r** is the boundary 
impedance, and Xm is the center of Dm, then the equations are derived, 
as a — > 0, for the effective field in the medium, consisting of many small 
particles, embedded in a given material (according to the recipe, derived in 
Section 2), under the following assumptions: 
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1) < K < 1, Ki = and (ESD holds, 

2) K> 1, Ki = i, and ([25]) holds. 

A remark is made about the case Ki = 1 and k = —1. In this case d = 0{a) 
is of the order of the size of a small particle. 

The results ar used for formulating a recipe for creating materials with a 
desired refraction coefficient. This recipe is similar to the one, given in [8], 
[5] in the case ki = 1/3 and k = 1. 
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